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(57)Abstract: 

PROBLEM TO BE SOLVED: To provide an intake air volume 
estimating device of an internal combustion engine for 
improving the calculation accuracy of the cylinder filling intake 
air volume by accurately calculating a spitting flow rate to an 
intake port from a cylinder by taking into consideration the 
pressure of the intake port. 

SOLUTION: When estimating an intake air flow rate mc filled in 
the cylinder by using a throttle valve model M10 t an intake pipe 
model M20, an intake valve model M30, a cylinder model M40, 
and a vent valve model M50 of modeling respective elements of 
an internal combustion engine intake system, when causing the 
spitting of making burnt gas in the cylinder flow backward to the 
intake port when opening/closing the intake valve, intake port 
pressure is calculated separately from intake pipe pressure by 
using an intake port model M25. Thus, since the accurate 
spitting gas volume is calculated by adding the influence of a 
local pressure increase in the intake port caused by spitting, 
the estimation accuracy of the cylinder filling intake air volume 
is improved. 
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* NOTICES * 

JPO and NCI PI are not responsible for any 
damages caused by the use of this translation. 

1 .This document has been translated by computer. So the translation may not reflect the original precisely. 
2.**** shows the word which can not be translated. 
3. In the drawings, any words are not translated. 



CLAIMS 



[Claim(s)] 
[Claim 1] 

In the inhalation air content presumption equipment of the internal combustion engine which computes the 
restoration inspired air volume in a cylinder using the inlet- valve passage inhalation-of-air flow rate which 
computed and computed the inlet- valve passage inhalation-of-air flow rate based on the inlet-valve model 
formula beforehand defined using the pressure-of-induction-pipe force and cylinder internal pressure, 
It has an inlet port pressure presumption means to compute a cylinder inlet port pressure based on the inlet 
port model formula beforehand defined using said inlet- valve passage inhalation-of-air flow rate, 
Inhalation air content presumption equipment of the internal combustion engine which presumes the 
restoration inspired air volume in a cylinder based on said inlet port pressure and cylinder internal pressure 
which were presumed. 
[Claim 2] 

Said inlet port pressure presumption means is inhalation air content presumption equipment of an internal 
combustion engine according to claim 1 which computes the value same when the blow return to which the 
residual gas in a cylinder passes an inlet valve, and flows backwards to an inlet port has not arisen as the 
pressure-of-induction-pipe force as an inlet port pressure, and computes the pressure which added a part for 
the pressure buildup depended for blowing over the pressure-of-induction-pipe force when blow return has 
arisen as an inlet port pressure. 
[Claim 3] 

the volume of the imagination which becomes settled according to engine operational status based on the 
relation defined beforehand when said blow return has produced said inlet port pressure presumption means 
— using — said ~ blowing — the inhalation air content presumption equipment of an internal combustion 
engine according to claim 2 which computes an inlet port pressure, assuming the pressure-buildup width of 
face in this space at the time of presupposing that the gas which flows backwards to an inlet port flowed in 
the space of said virtual volume to be a part for said pressure buildup which depends for blowing. 
[Claim 4] 

said virtual volume — further — said — blowing — the inhalation air content presumption equipment of an 
internal combustion engine according to claim 3 which changes based on the relation beforehand defined 
according to the amount of the gas which flows in this space at the time of presupposing that the gas which 
flows backwards to an inlet port flowed in the space of said virtual volume, and the pressure-buildup width 
of face in this space. 
[Claim 5] 

Said inlet port pressure presumption means is inhalation air content presumption equipment of an internal 
combustion engine according to claim 2 which computes a part for said pressure buildup depended for 
blowing based on the pressure loss in the inlet valve computed based on said inlet- valve passage inhalation- 
of-air flow rate when said blow return has arisen. 
[Claim 6] 

Said inlet port pressure presumption means is inhalation air content presumption equipment of an internal 
combustion engine according to claim 2 which computes a part for said inlet port pressure buildup depended 
for blowing based on the gas flow rate which flows into an inlet pipe, and said inlet- valve passage 
inhalation-of-air flow rate from the inlet port computed based on an inlet port pressure and the pressure-of- 
induction-pipe force when said blow return has arisen. 
[Claim 7] 

Said inlet port pressure presumption means is inhalation air content presumption equipment of an internal 
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combustion engine according to claim 2 which has the amplitude and period which are defined based on the 
relation which defined beforehand a part for said inlet port pressure buildup after inlet- valve valve opening 
according to engine operational status when said blow return has arisen, and is computed as a sine function 
which makes the starting point the time of inlet- valve valve-opening initiation. 



[Translation done.] 
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DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[Field of the Invention] 

This invention relates to the inhalation air content presumption equipment which computes an internal 
combustion engine's inhalation air content in a detail by the formula which modeled an internal combustion 
engine's induction system about an internal combustion engine's inhalation air content presumption 
equipment. 
[0002] 

[Description of the Prior Art] 

While dividing an internal combustion engine's induction system into elements, such as a throttle valve, an 
inlet pipe, an inlet valve, a cylinder, and an exhaust valve, modeling each element and expressing with a 
formula, the inhalation air content presumption equipment of the internal combustion engine using the so- 
called air model which calculates an engine's inhalation air content (restoration inspired air volume in a 
cylinder) by count is known by connecting between [ each ] models using a pressure, temperature, a flow 
rate, etc. 
[0003] 

With the inhalation air content presumption equipment using such an air model, it usually becomes possible 
to compute the restoration inspired air volume in a cylinder only by only whenever [ engine rotational 
frequency and throttle valve-opening ] in addition to atmospheric pressure and atmospheric temperature. 
Since the restoration inspired air volume in a cylinder is computable with the inhalation air content 
presumption equipment using such an air model with the count which used the model type, it is possible to 
compute an inhalation air content to responsibility fitness at the time of transient operation with the large 
change rate of whenever [ throttle valve-opening ] etc. 
[0004] 

As an example of the control unit of the internal combustion engine which performs this kind of inhalation 
air content presumption, there are some which were indicated by the patent reference 1 , for example. 
The equipment of the patent reference 1 calculates the pressure of the inhalation-of-air style in each model, 
temperature, and a flow rate using physical laws, such as the law of conservation of energy, a law of mass 
conservation, and a law of conservation of momentum, while it divides an engine's inhalation-of-air system 
into each element of a throttle valve, the inlet pipe containing a surge tank, and an inlet valve and expresses 
these elements with a simulation model. The restoration inspired air volume in a cylinder is computed based 
on the inhalation-of-air flow rate which passes the inlet valve computed by the above. 
[0005] 

With the equipment of the patent reference 1 , the restoration inspired air volume in a cylinder (namely, 
inspired air volume in a cylinder at the time of inlet-valve clausilium) is calculated as a value proportional to 
the pressure of the inlet pipe at the time of inlet- valve clausilium (surge tank), assuming that cylinder 
internal pressure is equal to the pressure-of-induction-pipe force (surge tank pressure) by the inhalation of 
air with which it filled up in the cylinder at the time of inlet-valve clausilium. 

However, with the equipment of the above-mentioned patent reference 1 , although the so-called blow return 
to which the gas in a cylinder flows backwards from an inlet valve to an inlet port with the pressure in a 
cylinder depending on an engine's service condition at the time of inlet-valve valve opening may arise in a 
actual engine, since this blow return is not taken into consideration, there is a problem to which the 
calculation precision of the restoration inspired air volume in a cylinder falls. 
[0006] 
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These people have already proposed an internal combustion engine's inhalation air content presumption 
equipment using the air model which took into consideration account the blow return of a top by the 
application for patent No. 316350 [ 2001 to ], in order to solve this problem. 

The above-mentioned application for patent No. 316350 [ 2001 to ], with public equipment, while dividing 
an engine's inhalation-of-air system into the element of a throttle valve, the inlet pipe containing a surge 
tank, and an inlet valve and modeling it, the air model which modeled the element of a cylinder and an 
exhaust valve further in addition to these is used. Although the throttle valve and the inlet pipe use the same 
model as the patent reference 1 , he is trying to compute the restoration inspired air volume in a cylinder in 
consideration of the blow return at the time of inlet- valve valve opening mentioned above with the 
equipment of this application for patent No. 3 16350 [ 2001 to ] by modeling an exhaust valve, a cylinder, 
and an inlet valve in addition to these, and computing the temperature of the inhalation of air in these 
elements, a pressure, and a flow rate using said each physical law. 
[0007] 

That is, with the equipment of the above-mentioned application for patent No. 316350 [ 2001 to ], the 
pressure in a cylinder is computed based on physical laws, such as relation between each model, a law of 
mass conservation, and the law of conservation of energy, and the inhalation-of-air flow rate which passes 
an inlet valve based on the pressure-of-induction-pipe force and cylinder internal pressure is computed. This 
flow rate turns into a forward flow rate which goes in a cylinder from an inlet pipe, when the pressure-of- 
induction-pipe force is higher than the pressure in a cylinder, but it becomes a negative value when the 
pressure in a cylinder is higher than the pressure-of-induction-pipe force (i.e., when blow return has arisen). 
The restoration inspired air volume in a cylinder is calculated by integrating this inlet- valve passage inspired 
air volume by the 1 cycle inlet-valve valve-opening period. 
[0008] 

[Patent reference 1 ] 

JP,2001-41095,A 

[0009] 

[Problem(s) to be Solved by the Invention] 

However, according to subsequent research, especially by having computed the inlet-valve passage 
inhalation-of-air flow rate only based on the pressure-of-induction-pipe force and cylinder internal pressure 
containing a surge tank like the equipment of the above-mentioned application for patent No. 3 16350 [ 2001 
to ], it has become clear that the calculation precision of the restoration inspired air volume in a cylinder 
falls greatly by a service condition which is large, i.e., conditions which are large. 
[0010] 

since the inlet-pipe volume which blows and contains a surge tank compared with a flow rate although the 
pressure-of-induction-pipe force computed will be blown and it will go up according to a flow rate with the 
equipment of the above-mentioned application for patent No. 316350 [ 2001 to ] if blow return arises is 
large — blowing — the lifting width of face of the pressure-of-induction-pipe force to produce becomes 
comparatively small, for this reason — blowing — a value with the comparatively large difference of the 
cylinder internal pressure at the time of generating, and the pressure-of-induction-pipe force — becoming — 
blowing — the amount of the gas which flows out of the inside of a cylinder size-comes comparatively, and 
comes. Thus, also when cylinder internal pressure falls on count when [ a lot of] it blows and a gas flows 
into an inlet pipe, and the inflow of the inhalation of air from an inlet port to into a cylinder is resumed, in 
order that the timing to which new mind flows in a cylinder may be overdue, the amount of the new mind 
with which it fills up in a cylinder decreases. 
[0011] 

however — a actual engine — the time of inlet- valve valve opening — blowing — the amount of the gas which 
flows backwards is smaller than the amount computed with the air model of the above-mentioned 
application for patent No. 316350 [ 2001 to ], and the restoration inspired air volume in a cylinder (new air 
volume) becomes larger than the value computed with the above-mentioned model. 
It is considered to be for a rapid local pressure buildup to arise into the inlet port part immediately after 
inlet-valve valve opening in a actual engine that blow in this way and a difference arises in an amount as a 
result of research. 
[0012] 

When blow return arises and especially bulb overlap is large, immediately after inlet-valve valve opening, 
the high speed almost near acoustic velocity blows from an inlet valve, and an air current flows into an inlet 
port. In the above-mentioned application for patent No. 316350 [ 2001 to ], it has calculated the inlet-valve 
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passage inhalation-of-air flow rate (blowing flow rate) by this high-speed flow having flowed into the inlet 
pipe (surge tank), and the pressure in a surge tank having assumed it to go up to homogeneity. For this 
reason, the pressure buildup of a surge tank becomes comparatively small, since the difference of a surge 
tank and cylinder internal pressure becomes comparatively large with as, it is blown, and a flow rate 
becomes comparatively large. 
[0013] 

However, in the inlet port where it blows immediately after inlet- valve valve opening actually, and this air 
current flows since the rate of flow of an air current is high-speed, a rapid and local pressure buildup 
happens and an inlet port pressure becomes quite higher than a surge tank pressure. For this reason, 
immediately after inlet-valve valve opening, the pressure differential of an inlet port and the inside of a 
cylinder contracts by lifting of inlet port internal pressure, and, thereby, the blow return by the inlet port 
from a cylinder comes to be controlled. 
[0014] 

That is, it will blow actually by the local pressure buildup in the inlet port immediately after inlet-valve 
valve opening, and a flow rate will decrease. 

For this reason, if it blows based on the pressure-of-induction-pipe force (surge tank pressure) and cylinder 
internal pressure and the flow rate is computed like the equipment of the above-mentioned application for 
patent No. 316350 [ 2001 to ], in a service condition with the large bulb overlap to which it blows especially 
and a flow rate becomes large, the problem to which the calculation precision of the restoration inspired air 
volume in a cylinder falls greatly will arise. 
[0015] 

In view of the above-mentioned problem, by taking the pressure of an inlet port into consideration, it blows, 
a flow rate is computed and this invention aims at offering the inhalation air content presumption equipment 
of the internal combustion engine in which the thing from a cylinder to an inlet port for which the 
calculation precision of the restoration inspired air volume in a cylinder is raised is possible to accuracy. 
[0016] 

[Means for Solving the Problem] 

According to invention according to claim 1 , based on the inlet- valve model formula beforehand defined 
using the pressure-of-induction-pipe force and cylinder internal pressure, an inlet- valve passage inhalation- 
of-air flow rate is computed. In the inhalation air content presumption equipment of the internal combustion 
engine which computes the restoration inspired air volume in a cylinder using the computed inlet- valve 
passage inhalation-of-air flow rate It has an inlet port pressure presumption means to compute a cylinder 
inlet port pressure based on the inlet port model formula beforehand defined using said inlet- valve passage 
inhalation-of-air flow rate. The inhalation air content presumption equipment of the internal combustion 
engine which presumes the restoration inspired air volume in a cylinder based on said inlet port pressure and 
cylinder internal pressure which were presumed is offered. 
[0017] 

That is, in invention of claim 1 , since an inlet port pressure presumption means to presume an inlet port 
pressure apart from the pressure-of-induction-pipe force using the inhalation-of-air flow rate which passes 
an inlet valve is established, for example, also when blow return arises, an inlet port pressure blows based 
on the inlet port pressure and cylinder internal pressure which it was presumed by accuracy and presumed, 
and a flow rate comes to be computed by accuracy. For this reason, it becomes possible to blow and to 
presume the restoration inspired air volume in a cylinder to accuracy also at the time of operation of bulb 
overlap size with many gaseous amounts. 
[0018] 

According to invention according to claim 2, said inlet port pressure presumption means When the blow 
return to which the residual gas in a cylinder passes an inlet valve, and flows backwards to an inlet port has 
not arisen The same value as the pressure-of-induction-pipe force is computed as an inlet port pressure, and 
when blow return has arisen, the inhalation air content presumption equipment of an internal combustion 
engine according to claim 1 which computes the pressure which added a part for the pressure buildup 
depended for blowing over the pressure-of-induction-pipe force as an inlet port pressure is offered. 
[0019] 

That is, in invention of claim 2, an inlet port pressure presumption means computes an inlet port pressure as 
a pressure higher than the pressure-of-induction-pipe force in consideration of the pressure buildup 
depended for blowing, when blow return has arisen. Thereby, also when blow return has arisen, the amount 
of return [ blow ] comes to be controlled compared with the case where it is the former assumed that the 
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pressure-of-induction-pipe force and an inlet port pressure are equal, and it becomes possible to presume the 

restoration inspired air volume in a cylinder to accuracy. 

[0020] 

Moreover, the rate of the inhalation-of-air style which flows in a cylinder through an inlet port in the usual 
inhalation-of-air style which blow return has not produced is blown, and is far smaller than an air current. 
For this reason, a big pressure differential is not produced between an inlet pipe (surge tank) and an inlet 
port. Therefore, when blow return has not arisen, it is not necessary to compute an inlet port pressure 
separately, and the restoration inspired air volume in a cylinder can be computed using the value same as an 
inlet port pressure as the pressure-of-induction-pipe force. 
[0021] 

According to invention according to claim 3, said inlet port pressure presumption means When said blow 
return has arisen, the volume of the imagination which becomes settled according to engine operational 
status based on the relation defined beforehand is used, said — blowing — the inhalation air content 
presumption equipment of an internal combustion engine according to claim 2 which computes an inlet port 
pressure, assuming the pressure-buildup width of face in this space at the time of presupposing that the gas 
which flows backwards to an inlet port flowed in the space of said virtual volume to be a part for said 
pressure buildup depended for blowing is offered. 
[0022] 

That is, in invention of claim 3, when it assumes that it is the closed volume section with the volume which 
becomes settled according to engine operational status about an inlet port, and blows, it blows over this 
closed volume section at the time of generating and a gas flows, it is assumed that the pressure buildup of an 
inlet port arises. Thereby, the pressure buildup of an inlet port can be computed according to the quantity of 
gas of return [ blow ], and it becomes possible to compute the restoration inspired air volume in a cylinder to 
accuracy. 
[0023] 

according to invention according to claim 4 — said virtual volume — further ~ said — blowing — the 
inhalation air content presumption equipment of an internal combustion engine according to claim 3 which 
changes based on the relation beforehand defined according to the amount of the gas which flows in this 
space at the time of presupposing that the gas which flows backwards to an inlet port flowed in the space of 
said virtual volume, and the pressure-buildup width of face in this space is offered. 
[0024] 

That is, in invention of claim 4, like claim 3, when blow return has arisen, it assumes that it is the volume 
section which closed the inlet port, and flows and blows over this volume section, and the pressure buildup 
of an inlet port is computed according to quantity of gas. However, in invention of claim 4, the volume of 
this closed volume section is not fixed, and it computes an inlet port pressure by assuming [ that it can 
expand and ] it according to the internal gas pressure force, namely, the virtual volume of an inlet port 
assuming it to change according to the pressure buildup of an inlet port. For example, the gas inside the 
volume section can assume expansion of the volume section to resist resistance of elastic resistance (spring), 
frictional resistance (dash pot), etc., and to expand the volume section. Thus, since it becomes possible to 
blow over accuracy further and to compute quantity of gas by modeling an inlet port as the volume section 
of the adjustable volume, it becomes possible more to compute the restoration inspired air volume in a 
cylinder to accuracy. 
[0025] 

According to invention according to claim 5, when said blow return has produced said inlet port pressure 
presumption means, the inhalation air content presumption equipment of an internal combustion engine 
according to claim 2 which computes a part for said pressure buildup depended for blowing based on the 
pressure loss in the inlet valve computed based on said inlet-valve passage inhalation-of-air flow rate is 
offered. 
[0026] 

That is, based on the inhalation-of-air flow rate (namely, back run flow rate depended for blowing) which 
blows by invention of claim 5 and sometimes passes an inlet valve, the pressure drawdown (pressure loss of 
the passage of imagination containing an inlet valve and the inlet port part near an inlet valve) in an inlet 
valve is computed, and a part for an inlet port pressure buildup is computed based on this pressure loss. 
Since it blows, it sometimes blows by this and an inlet port pressure is computed according to a flow rate, it 
blows, a flow rate is computed by accuracy and the calculation precision of the restoration inspired air 
volume in a cylinder improves. 
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[0027] 

According to invention according to claim 6, when said blow return has produced said inlet port pressure 
presumption means, the inhalation air content presumption equipment of an internal combustion engine 
according to claim 2 which computes a part for said inlet port pressure buildup depended for blowing based 
on the gas flow rate which flows into an inlet pipe, and said inlet-valve passage inhalation-of-air flow rate 
from the inlet port computed based on an inlet port pressure and the pressure-of-induction-pipe force is 
offered. 
[0028] 

That is, in invention of claim 6, an inlet port pressure is computed based on a law of mass conservation 
using the quantity of gas (that is, an inlet valve is passed and it flows backwards to an inlet port blowing 
flow rate) which flows an inlet port pressure when blow return has arisen into an inlet port, and the quantity 
of gas (quantity of gas which flows backwards from an inlet port to a surge tank) which flows out of an inlet 
port. The quantity of gas which flows into an inlet port, and the quantity of gas which flows out of an inlet 
port can be calculated as a gas flow rate which passes drawing based on an inlet port pressure, the cylinder 
pressure force, and an inlet port pressure and the pressure-of-induction-pipe force, respectively. Thereby, 
since an inlet port pressure is computed by accuracy, the calculation precision of the restoration inspired air 
volume in a cylinder improves. 
[0029] 

According to invention according to claim 7, when said blow return has produced said inlet port pressure 
presumption means, it has the amplitude and period which are defined based on the relation which defined 
beforehand a part for said inlet port pressure buildup after inlet- valve valve opening according to engine 
operational status, and the inhalation air content presumption equipment of an internal combustion engine 
according to claim 2 computed as a sine function which makes the starting point the time of inlet-valve 
valve-opening initiation is offered. 
[0030] 

That is, in invention of claim 7, the temporal response of the inlet port pressure buildup after inlet-valve 
valve opening when blow return has arisen is assumed to be a sine function, for example, it asks for the 
amplitude of this function, and a period beforehand based on the experimental value as a function of engine 
operational status (for example, an engine engine speed, a load, the amount of bulb overlap, etc.). By 
computing an inlet port pressure using the function based on such an actual measurement, since [ near a 
more nearly actual value / being exact ] it blows and a flow rate is computed, the calculation precision of the 
restoration inspired air volume in a cylinder improves. 
[0031] 

[Embodiment of the Invention] 

Hereafter, the operation gestalt of this invention is explained using an accompanying drawing. 
Drawing 1 is drawing showing the outline configuration of the operation gestalt at the time of applying the 
inhalation air content presumption equipment of this invention to a jump-spark-ignition type engine. 
[0032] 

In drawing 1 , a reference mark 10 shows a jump-spark-ignition type internal combustion engine body. With 
this operation gestalt, although the 4-cy Under gasoline engine is used as an internal combustion engine 1 0, 
this invention is not necessarily limited to a 4-cylinder gasoline engine, and can be applied also to the 
internal combustion engine of other formats. 
[0033] 

An engine 10 has a cylinder 21, a piston 22, and a crankshaft 24, and the combustion chamber 25 is formed 
in the piston 22 upper part in a cylinder 21. 

Moreover, the phase angle of the inhalation-of-air cam shaft which drives the inlet valve 32 which opens 
and closes the inlet port 31 which is open for free passage to a combustion chamber 25, and an inlet port 31 
to the cylinder head of the cylinder 21 upper part, and an inlet valve 32, and the amount of valve lifts of an 
inlet valve 32 are prepared in the inlet-valve control unit 33 which can be changed continuously. It is the 
actuator of the inlet- valve control device 33 which 33a shows to drawing 1 . Moreover, the fuel injection 
valve 39 which injects a fuel is formed in the inlet port in the inlet port 31 . 
[0034] 

Moreover, the exhaust air cam shaft 36 which carries out closing motion actuation of the exhaust valve 35 
which opens and closes an inlet port 31, the exhaust port 34 which is open for free passage to a combustion 
chamber 25 similarly, and an exhaust port, and the exhaust valve 35 is formed in the cylinder head of the 
cylinder 21 upper part. 
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Furthermore, the firing circuit (ignitor) 3 8 which faces the cylinder head of the cylinder 2 1 upper part in a 
combustion chamber 25, and supplies high tension to an ignition plug 37 and an ignition plug is formed. 
[0035] 

An engine's 10 induction system 40 is equipped with the surge tank 44 formed in the inlet pipe 41 connected 
to the inlet port 3 1 of each cylinder of an engine through the inlet manifold, and the inlet pipe 41, the air 
cleaner 42 formed in inlet-pipe 41 inlet port, and the throttle valve 43 arranged at the surge tank 44 upstream 
of an inlet pipe 41. 
[0036] 

With this operation gestalt, the throttle valve 43 is equipped with actuator 43 a of proper formats, such as a 
stepper motor, and let it be the so-called electronic throttle valve which can take the independent opening 
with accelerator pedal 82 actuation of an operator. 
[0037] 

Moreover, it is inhaled by the inlet pipe 41 of the throttle-valve 43 upstream at an inlet pipe, and the air flow 
meter 61 which detects the inhalation air content which passes a throttle valve 43 and flows is formed. With 
this operation gestalt, the air flow meter 6 1 is used as the hotwire flow meter, and the intake-air-temperature 
sensor 62 which detects an intake-air temperature, and the atmospheric pressure sensor 63 which detects 
atmospheric pressure are formed in about 61 air flow meter. 
[0038] 

Moreover, an engine's 10 exhaust system 50 is equipped with the exhaust pipe 52 connected to the exhaust 
port 34 of each cylinder through the exhaust manifold 51 . What 53 shows to drawing 1 is the catalytic 
converter which contained the exhaust air clarification catalyst. 

What 80 shows to drawing 1 is an electronic control unit (ECU) which controls an engine 10. ECU80 
performs exact inhalation air content presumption which was the microcomputer of a well-known 
configuration of having connected the backup RAM84 grade which can hold the content of storage by the 
bidirectional bus even if CPU81, ROM82, RAM83, the input/output interface 85, and the power source were 
disconnected, and basic control of engines, such as ignition timing and fuel-injection control, was 
performed, and also took blow return into consideration based on the air model later mentioned with this 
operation gestalt. 
[0039] 

The signal corresponding to an intake-air temperature (atmospheric temperature) and atmospheric pressure 
in the signal corresponding to the inhalation air content (mass flow) which flows into a throttle valve 43 
from the heat ray type air flow meter 61 is inputted into the input/output interface 85 of ECU80 from the 
intake-air-temperature sensor 62 and the atmospheric pressure sensor 63 again, respectively for these 
control. The output of the intake-air-temperature sensor 62 and the atmospheric pressure sensor 63 is used in 
order to amend the inhalation air content detected with the air flow meter 61 according to atmospheric 
temperature and atmospheric pressure. 

In addition, although not explained to a detail, with this operation gestalt, the inhalation air content detected 
with the air flow meter 61 is used for amending the inhalation air content computed using the air model 
mentioned later. 
[0040] 

The voltage signal corresponding to an operator's amount of accelerator pedal treading in (accelerator 
opening) is inputted, respectively from the accelerator opening sensor 71 by which the voltage signal 
corresponding to throttle- valve 43 opening was formed in about 82 accelerator pedal of an engine from the 
sensor 64 whenever [ throttle valve-opening / which was prepared in about 43 throttle valve ], and also the 
pressure-of-induction-pipe force is inputted into the input/output interface 85 of ECU80 from the intake- 
pressure sensor 65 arranged at the surge tank 44. 
[0041] 

furthermore — this operation gestalt — 02 sensor 70 is formed in the catalytic-converter 53 upstream of an 
exhaust pipe 52, the circulating-water-temperature sensor 69 is formed in an engine's 10 cooling water path, 
respectively, and the inlet-valve lift sensor 67 has inputted the amount of lifts, exhaust air oxygen density, 
and engine circulating water temperature of an inlet valve 32 into the input/output interface 85 of ECU80 at 
the inlet- valve control device 33, respectively. 
[0042] 

The crank angle sensor arranged near the crankshaft shows that a sign 66 shows to drawing 1 by the cam 
angle sensor arranged near the inhalation-of-air cam shaft, and 68. 

The cam angle sensor 66 outputs a cam position signal, whenever 180 degrees of crankshafts 24 rotate, and 
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the crank angle sensor 68 outputs a criteria crank location pulse signal with pulse width wider than the 
above-mentioned crank angle-of-rotation pulse signal, whenever a crankshaft 24 takes a criteria revolution 
location, while outputting a crank angle-of-rotation pulse signal, whenever a crankshaft 24 carries out a 
fixed include-angle (for example, 10 degrees) revolution (whenever the 1st cylinder arrives at a top dead 
center). 
[0043] 

These signals are inputted into the interface 85 of ECU80. ECU80 computes the revolution location (crank 
angle) of a crankshaft from the number of the crank angle-of-rotation pulse signals after a criteria crank 
location pulse signal input while computing engine rotational speed (rotational frequency) from spacing of 
for example, a crank angle-of-rotation pulse signal. Moreover, the cam position signal inputted from the cam 
angle sensor 66 is used for computing the valve timing of an inlet valve with a crank angle. 
[0044] 

Moreover, with this operation gestalt, the exhaust gas pressure Pe of an exhaust port, the exhaust-gas- 
temperature sensor which detects temperature Te, and the exhaust-pressure sensor are formed in the exhaust 
port of a cylinder, respectively (not shown), and the signal corresponding to an exhaust-gas temperature Pe 
and exhaust gas pressure Te is inputted into the interface 85 of ECU80, respectively. 
[0045] 

With this operation gestalt, ECU80 computes the fuel oil consumption fc in the cylinder as fc=KIxMc using 
the inspired air volume Mc with which an engine's 10 cylinder 20 is filled up. KJ is a conversion factor 
which becomes settled according to an engine air- fuel ratio. 
[0046] 

Moreover, ECU80 computes the inspired air volume (restoration inspired air volume in a cylinder) Mc with 
which a cylinder 20 is filled up using the air model formula which modeled an engine's 10 induction system 
40. 

[0047] 

With this operation gestalt, the applicant for this patent uses fundamentally the restoration inspired-air- 
volume calculation approach in a cylinder based on the air model proposed by the application for patent No. 
316350 [ 2001 to ]. Therefore, in order to explain the restoration inspired-air- volume Mc calculation 
approach in a cylinder of this operation gestalt, the restoration inspired-air- volume Mc calculation approach 
in a cylinder in an application for patent No. 316350 [ 2001 to ] is explained briefly first. 
[0048] 

An engine's 10 induction system 40 is disassembled into a throttle valve, the inlet pipe containing a surge 
tank, an inlet valve, a cylinder, and five elements of an exhaust valve, and he expresses each element with a 
simulation model, and is trying to compute the inspired air volume (restoration inspired air volume in a 
cylinder) Mc with which it fills up in a cylinder eventually with this operation gestalt by mathematizing the 
pressure of the inhalation of air which flows the inside of an element, temperature, and a flow rate. 
[0049] 

Hereafter, each model is explained. 
(1) Throttle- valve model 

A throttle- valve model is given as flow of the gas which passes drawing (throttle valve-opening regio oralis) 

as fundamentally shown in drawing 2 . 

[0050] 

That is, as shown in drawing 2 , when a throttle valve 43 is considered as drawing, the gas flow rate mt 
which passes a throttle valve 43 as known well becomes settled with the gas pressure force before and 
behind a throttle valve, temperature, and throttle-valve 43 opening area, and is expressed with the throttle- 
valve model formula of the following several 1 or several 2. 
[0051] 
[Equation 1] 



' V2KRTa VU + V I Pa k + U 
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[0052] 
[Equation 2] 



rat = Ct (9t ) At (0t) PmJ K + 1 Jf-M ' - [ — — 1 

v ' v ' V2KRTm vU + lJ \Ym K + lJ 



[0053] 

Here, several 1 is carrying out the table of the throttle- valve passage air flow rate when the pressure of the 
throttle- valve 43 upstream is lower than the pressure of the downstream and, as for several 2, the back run 
has produced the throttle-valve passage inhalation-of-air flow rate of the usual direction of inhalation-of-air 
flow (forward feed) when the pressure of the throttle-valve 43 upstream (atmospheric-air side) is higher than 
the pressure of the downstream (inlet-pipe side), respectively. 
[0054] 

In several 1 and several 2, Ct (thetat) is a flow coefficient, At (thetat) is a throttle valve opening area, and Ct 
(thetat) and At (thetat) are the functions of thetat whenever [ throttle valve-opening ], respectively. Pa is the 
pressure of the throttle- valve upstream, and Pm expresses atmospheric pressure with this operation gestalt, 
and it expresses the pressure-of-induction-pipe force with this operation gestalt by the pressure of the 
throttle- valve downstream. Moreover, Ta and Tm are the intake-air temperatures of the throttle- valve 
upstream and the downstream, respectively, atmospheric temperature and whenever [ inlet-pipe internal 
temperature ], R expresses a gas constant and kappa expresses the ratio of specific beat of inhalation of air 
with this operation gestalt, respectively. In addition, with this operation gestalt, R and kappa are treated as 
constant value. 
[0055] 

In this operation gestalt, a flow coefficient Ct (thetat) and the throttle valve opening area At (thetat) are 
beforehand computed by every theta whenever [ throttle valve-opening ], and are stored in ROM82 of 
ECU80 in the form of a numerical map. 

As shown in several 1 and several 2, the throttle-valve model of this operation gestalt is a model formula 
which considers [ whenever / atmospheric pressure Pa and pressure-of-induction-pipe force Pm, 
atmospheric temperature Ta or inlet-pipe internal temperature ] the throttle-valve passage inhalation-of-air 
flow rate mt as an output by considering thetat as an input whenever [ Tm and throttle valve-opening ]. 
[0056] 

(2) Inlet-pipe model 

Drawing 3 shows an inlet-pipe model. An inlet-pipe model computes the pressure-of-induction-pipe force 
Pm and temperature Tm using a law of mass conservation and the law of conservation of energy based on 
the inhalation-of-air flow rate which flows in the space of the fixed volume, and the inhalation-of-air flow 
rate flowing out. 
[0057] 

Since the inhalation-of-air style which flows into the volume section Vm is a throttle-valve passage 
inhalation-of-air style as shown in drawing 3 when the volume of the inlet pipe 41 containing the surge tank 
44 from a throttle valve 43 to an inlet valve 32 is set to Vm, an inflow inhalation-of-air flow rate serves as 
mt, and an intake-air temperature serves as Ta. 
[0058] 

Moreover, the inhalation-of-air flow rate which flows out of the volume section Vm is the inhalation-of-air 
flow rate mc with which passes an inlet valve and a cylinder is filled up, and the temperature of this 
inhalation-of-air style turns into the inlet-pipe temperature Tm. 
[0059] 

Therefore, the four following are drawn for the three following from the law of conservation of energy from 
a law of mass conservation, respectively. 
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[0060] 
[Equation 3] 



d (Pm / Tm ) / dt = (R / Vm) • (mt - mc) 



[0061] 
[Equation 4] 



dPm/dt = K-(R/Vm) (mt Ta-mc-Tm) 



[0062] 

With this operation gestalt, in order to perform data processing in ECU80, the following several 5 and 
several 6 formulas which integrated with several 3 and several 4 and were discretized are used. 
[0063] 
[Equation 5] 



^(k)= ~(k- 1)+ At-5-(mt(k - l)-mcAVE(k - 1)) 
1 m 1 m Vm 



[0064] 
[Equation 6] 



Pm (k)= Pm (k - l)+AtK-^-(mt(k -l)Ta- mcAVE (k - l)Tm (k - 1)) 

Vm 



[0065] 

Here, Subscript k expresses the value at the time of count activation, k-1 expresses the value at the time of 
count activation last time this time, and deltat expresses between the difference minutes at the time of being 
discretized (for example, mesh for about 8ms). Moreover, mcAVE in several 5 and several 6 formulas: It is 
the time amount average crossed to a part for the 1 cycle intake stroke of the inhalation-of-air flow rate mc 
in a cylinder mentioned later. 
[0066] 

The inlet-pipe temperature Tm is computed as Tm=Pm/(Pm/Tm) using Pm/Tm calculated from several 5, 
and Pm calculated from several 6. 
(3) Inlet-valve model 

An inlet-valve model is a model showing the inhalation-of-air style which passes an inlet valve 32, as shown 
in drawing 4 , and it is fundamentally expressed as a throttle- valve model as an inhalation-of-air style which 
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passes drawing (inlet-valve opening area) similarly. 
[0067] 

Therefore, an inlet-valve model is expressed with the several 7 and several 8 following formula like the 
above-mentioned throttle- valve model (several 1, several 2). 
[0068] 
[Equation 7] 



mc = Civ(L)Aiv(L)PmJ— ^±J_ 

2)cRTm 




[0069] 
[Equation 8] 



mc = Civ(L)Aiv(L)Pc/^±l- 1 -J^Y-f^-_iJ 

' v ' V2kRTc vU + lJ I, Pc K + lJ 



[0070] 

Here, the pressure-of-induction-pipe force and temperature, and Pc and Tc of the inhalation-of-air flow rate 
to which mc passes an inlet valve, i.e., the inspired air volume with which it fills up in unit time amount in a 
cylinder, and Pm and Tm are cylinder internal pressure and temperature. Moreover, Civ (L) is a flow 
coefficient, and the function of the amount L of lifts of an inlet valve and Aiv (L) are the opening area of an 
inlet valve 32, and it serves as a function of the amount L of lifts similarly. 
[0071] 

In this operation gestalt, the product value Civ (L) and the value of Aiv (L) are calculated to L value each, 
and are stored in ROM82 of ECU80 in the form of a numerical table. 

Like the throttle- valve model (several 1 , several 2), several 7 has the pressure-of-induction-pipe force Pm 
larger than the cylinder internal pressure Pc, and is carrying out the table of the case where several 8 Pm is 
smaller than Pc, and the inhalation of air out of a cylinder is flowing backwards the case where inhalation of 
air is flowing in a cylinder, to the inlet pipe, respectively. 
[0072] 

As shown in several 7 and several 8, the inlet- valve model of this operation gestalt outputs [ whenever / 
pressure-of-induction-pipe force Pm, cylinder internal pressure Pc, and inlet-pipe internal temperature ] the 
inlet- valve passage inhalation-of-air flow rate mc whenever [ Tm or cylinder internal temperature ] by 
considering Tc and the amount L of inlet-valve lifts as an input. 

The five above-mentioned inlet-pipe model formulas and the average inhalation-of-air flow rate mcAVE 
used by several 6 are the values which **(ed) the value Mc with which it integrated over the whole inlet- 
valve passage air flow rate mc intake stroke computed by several 7 and several 8 by the time amount T 
(180CA) which is equivalent to 180 degrees of crank angles in **. 
[0073] 

Namely, mcAVE=Mc/T (180CA) 
(4) Cylinder model 

Although the cylinder model of this operation gestalt is expressed as a model with which the inhalation of 
air of a flow rate mc and temperature Tm flows into the closed space (Vc is the function of a crank angle) of 
the volume Vc (adjustable) per unit time amount from an inlet valve fundamentally as shown in drawing 5 It 
is given as a model with which exhaust air of a flow rate me and temperature Te flows during inlet-valve 

http://www4.ipdl .ncipi . go j p/cgi-bin/tran_web_cgi_ejj e 6/8/2006 



JP,2004-278359,A [DETAILED DESCRIPTION] 



Page 1 1 of 29 



valve opening in consideration of the flow of the exhaust air which passes an exhaust valve from an exhaust 
port and flows backwards in a cylinder in addition to the inhalation-of-air style from an inlet valve. 
[0074] 

That is, in the cylinder model of this operation gestalt, mc, me, Tm, and Te (and crank angle) are considered 
as an input, and the total amount Mc of the inhalation of air with which it filled up in the cylinder into the 
cylinder internal pressure Pc, temperature Tc, and an intake stroke is outputted. 
[0075] 

A cylinder model is fundamentally expressed with the nine following based on the law of conservation of 

energy. 

[0076] 

[Equation 9] 



dPc KRTm tcRTe kPc dVc „ _Qw 

= # mc+ •me • mk - 1) 

dt Vc Vc Vc dt Vc 



[0077] 

Here, Qw is a heating value given to inhalation of air from a cylinder wall, an inlet port wall surface, etc. 
[0078] 

Moreover, Tc is expressed with the ten following from a gaseous equation of state whenever [ cylinder 

internal temperature ]. 

[0079] 

[Equation 10] 



Tc= (PcVc)/(MclR) 



[0080] 

Here, Mcl is the restoration inspired air volume from the time of inlet- valve valve opening, and is given as 
an integral value of the inflow gas flow rate (mc+me) from the inlet valve and exhaust valve from the time 
of inlet-valve valve opening to current (at the time of Tc count). 
[0081] 

By the form of this as, since nine above cannot be calculated by ECU80, it is changed to a discretization 

format and used in the following several 1 1 forms. 

[0082] 

[Equation 1 1 ] 
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[0083] 

In addition, deltat in several 1 1 is between the difference minutes at the time of several 9 being discretized 
(in a cylinder model, compared with the case of an inlet-pipe model, deltat is set up short, for example, is set 
as the mesh for about 1ms). 
[0084] 

Moreover, the exhaust valve passage flow rate me is calculated with the exhaust valve model explained 
below. 

(5) Exhaust valve model ( drawing 7 , M50) 

An exhaust valve model is a model which outputs the exhaust valve passage flow rate me like an inlet-valve 
model by considering the pressure, the temperature, Pe and Te and the cylinder pressure force, the 
temperature Pc and Tc, and the crank angle of an exhaust port as an input. 
[0085] 
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Drawing 6 shows an exhaust valve model. As shown in drawing 6 , the exhaust valve model is 
fundamentally the same as an inlet-valve model, and an exhaust valve passage flow rate is given by several 
1 2 or several 1 3 . 
[0086] 

[Equation 12] 



me = Cev(L)- Aev(L)*Pe 



CUT f -JS-Y-fJ*- — 1 

V 2xRTe y V k + 1 J \Pe K + lJ 



[0087] 

[Equation 13] 



Several 1 2 shows the case where exhaust air flows from several 1 3 cylinder to an exhaust port, respectively, 
when exhaust air flows backwards in a cylinder from an exhaust port (when Pc is higher than Pe). (namely, 
when the cylinder internal pressure Pc is lower than the exhaust port pressure Pe) In addition, Cev (L) and 
Aev (L) are the flow coefficients and opening area of an exhaust valve 35, respectively, and both serve as a 
function of the amount L of lifts of an exhaust valve. 
[0088] 

In addition, although the exhaust port pressure Pe and an exhaust-gas temperature Te are detected with this 
operation gestalt by the exhaust-pressure sensor and exhaust-gas-temperature sensor (not shown) which 
were formed in the exhaust port, since it is comparatively small, the pressure loss of an exhaust air system 
can also be calculated by assuming Pe to be equal to atmospheric pressure in approximation instead of an 
exhaust-pressure sensor detecting an exhaust pressure. 

Moreover, an exhaust-gas temperature Te considers abbreviation etc. as whenever [ cylinder bashful 
temperature ], and is hurt, and you may make it substitute Tc whenever [ cylinder internal temperature ] 
similarly, without forming an exhaust-gas-temperature sensor also about an exhaust-gas temperature Te. 
[0089] 

Drawing 7 is the block diagram showing the whole air model which consists of the throttle-valve model 
explained above, an inlet-pipe model, an inlet- valve model, a cylinder model, and an exhaust valve model. 
The model formula of each model element, an input and an output, the relation between each model 
element, the calculation process of the restoration inspired air volume Mc in a cylinder, etc. are again 
summarized using drawing 6 . 
[0090] 

(A) Throttle- valve model ( drawing 7 , M10) 

( 1 ) Formula 

[0091] 

[Equation 14] 
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mt-Ct(et)At(9t)Pa K+l " K 1 ^ Pm 




2KRTa VVK + ly ^Pa k + 1J 



[0092] 

[Equation 1 5] 



mt = Ct(0t)At(Gt)Pm K + 1 



2icRTm 




[0093] 

(2) Input 

- Throttle valve-opening whenever thetat 

- Atmospheric pressure Pa (atmospheric pressure sensor 63 output) 

- Atmospheric temperature Ta (intake-air- temperature sensor 62 output) 

- Pressure-of-induction-pipe force Pm 

- Inlet-pipe temperature Tm (in the case of several 1 5 (at the time of a back run)) 

(3) Output 

- Throttle-valve passage inhalation-of-air flow rate mt 
(B) Inlet-pipe model ( drawing 7 , M20) 

(1) Formula 
[0094] 

[Equation 16] 



^00= ^(k- 1) + At^-(mt(k - 1)- mcAVE(k - 1)) 
Tm Tm Vm 



[0095] 

[Equation 17] 



Pm (k)= Pm (k - 1)+ AtK-5-(mt (k -l)Ta- mcAVE(k - l)Tm (k - 1)) 

Vm 



[0096] 

[Equation 18] 
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Tm(k)= Pm < k ) Pni 



[0097] 

(2) Input 

- Throttle- valve passage inhalation-of-air flow rate mt (k-1) 

- Inlet- valve passage inhalation-of-air flow rate mcAVE (average) (k-1) 

- Pm (k-1), Tm (k-1) (k (Pm/Tm)-1), 

(3) Output 

- Pressure-of-induction-pipe force Pm (k) 

- Inlet-pipe temperature Tm (k) 

(C) Inlet- valve model ( drawing 7 , M30) 

(1) Formula 

[0098] 

[Equation 19] 



mc = Civ (L) Aiv (L) Pm 




(— -— l 

\,Pm K + lJ 



[0099] 

[Equation 20] 



mc = Civ(L)Aiv(L)Pc 



K + l 

2kRTc 



f K ] 


2 


( Pm 






Ipc 




[0100] 

(2) Input 

- Pressure-of-induction-pipe force Pm 

- Inlet-pipe temperature Tm 

- Cylinder pressure force Pc 

- It is Tc (in the case of several 20 (at the time of a back run)) whenever [ cylinder internal temperature ]. 

- The amount L of valve lifts (from a crank angle to calculation) 

(3) Output 

Inlet- valve passage inhalation-of-air flow rate mc (instantaneous value) 

- Inlet-valve passage inhalation-of-air flow rate mcAVE (average) 
(D) Cylinder model ( drawing 7 , M40) 

(1) Formula 
[0101] 

[Equation 2 1 ] 
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[0102] 

[Equation 22] 



Tc= (Pc Vc)/(Mcl R) 



[0103] 
(2) Input 
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- Cylinder capacity Vc (from a crank angle to calculation) 

- Inlet-valve passage inspired air volume mc (instantaneous value) 

- It is Tm whenever [ inlet-pipe internal temperature ]. 

- Exhaust valve passage inspired air volume me 

- Exhaust-gas temperature Te 
(3) Output 

- Cylinder pressure force Pc 

- It is Tc whenever [ cylinder internal temperature ]. 

- Restoration inspired air volume Mc in a cylinder 
(E) Exhaust valve model ( drawing 7 , M50) 

(1) Formula 
[0104] 

[Equation 23] 

me = Cev(L)-Aev(L)-Pe.f-^- Jf-^l -f— — 1 



[0105] 

[Equation 24] 



me = Cev(L)- Aev(L)«Pc 



/^±r|r_«_y_ri»__L.Y 

V2kRTc ^K + lJ \^Pc K + lJ 



[0106] 

(2) Input 

- Exhaust pipe pressure Pe 

- Exhaust pipe temperature Te 

- Cylinder pressure force Pc 

- It is Tc (in the case of several 24) whenever [ cylinder internal temperature ]. 

- The amount L of valve lifts (from a crank angle to calculation) 

(3) Output 

Exhaust valve passage inhalation-of-air flow rate me (instantaneous value) 

As mentioned above, the input of models other than a throttle-valve model and an exhaust valve model, i.e., 
an inlet-pipe model, an inlet-valve model, and a cylinder model can use the output value of other models 
except [ all ] a crank angle. Moreover, a crank angle is computable as a function of the engine rotational 
frequency NE. 
[0107] 

Furthermore, the input of a throttle- valve model can also use the output value of other models except for 
throttle valve-opening whenever theta (t), an atmospheric pressure Pa, and atmospheric temperature Ta. And 
since atmospheric pressure Pa and atmospheric temperature Ta do not change, they may be considered to be 
abbreviation constant value for a short time. 

Therefore, in the air model shown in drawing 7 , if throttle valve-opening whenever theta (t), the engine 
rotational frequency NE and exhaust gas pressure Pe, and temperature Te become settled fundamentally, the 
restoration inspired air volume mc in a cylinder is computable. 
[0108] 
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concrete — for example, the time of engine start up — as initial value — Pm and Pc — atmospheric pressure 
Pa — Tm and Tc — atmospheric temperature Ta — moreover, mc is set as flow rates mt and 0, respectively, 
and above-mentioned several 14 to 24 is calculated using throttle valve-opening whenever theta (t) and a 
crank angle (Pe and Te use an actual measurement). 

And count of several 14 to several 24 next time is performed using values, such as Pm and PcTm which 
were obtained as a result, and Tc, mt, mc, me. And this count is repeated using the value of theta at that 
event (t), and NE, and values, such as Pm, tm, mt, mc, me, etc. which were calculated by the last count, at 
predetermined spacing (for example, about 8ms) from the time of engine start up. Thereby, it converges on a 
actual value in a short time after count initiation, and Pm, tm, mt, mc, etc. must come to be in agreement 
with the actual restoration inspired air volume in a cylinder. 
[0109] 

In the above-mentioned conventional air model, in the inlet- valve model ( drawing 4 ), since inspired-air- 
volume presumption also in consideration of the case (several 20) where the back run to an inlet port arises 
from a cylinder is performed, for example, also when a lot of [ the amount of bulb overlap is large and ] 
blow return arises, the restoration inspired air volume in a cylinder should be presumed by accuracy. 
However, actually, as mentioned above, since the conventional air model is blown and is not taking into 
consideration the local pressure buildup in the inlet port at the time, the problem to which it blows over, and 
the presumed precision of the restoration inspired air volume in a cylinder falls especially when an amount 
is large produces it. 
[0110] 

In this invention, the presumed precision of the restoration inspired air volume in a cylinder is raised by 
dividing and modeling an inlet pipe and an inlet port to the above-mentioned conventional air model, 
blowing using an inlet port model, and computing the inlet port pressure at the time of generating as a value 
with the separate pressure-of-induction-pipe force. 

Drawing 8 is the same block diagram as drawing 7 explaining the whole air model configuration of this 

operation gestalt. 

[0111] 

As shown in drawing 8 , with this operation gestalt, the point that the inlet port model M25 is formed 
between the inlet-pipe model M20 and the inlet- valve model M30 is different. 

Although later mentioned about the detail of the inlet port model M25, he computes the inlet port pressure 
Pp using this inlet port model M25, and is trying to compute the inlet- valve passage flow rate mc using 
several 25 shown below instead of several 20, and several 26 with this operation gestalt. Here, several 25 is 
the case where blow return has produced several 26 in the inlet port, when the blow return by the inlet port 
has not arisen (namely, when it is Pp<Pc). (namely, when it is Pp>=Pc) 
[0112] 

[Equation 25] 



mc = Civ(L)Aiv(L)Pp K + 1 



2KRTm \H k + 1 ) 





[0113] 

[Equation 26] 
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[0114] 

When the blow return from a cylinder has arisen, with the inlet port model M25 of this operation gestalt, the 
inlet port pressure Pp is computed as a pressure higher than Pm separately from the pressure-of-induction- 
pipe force Pm here, so that it may mention later, but when blow return has not arisen, a value equal to the 
pressure-of-induction-pipe force Pm as an inlet port pressure Pp is used, and it does not perform computing 
the inlet port pressure Pp separately from the pressure-of-induction-pipe force Pp. 
[0115] 

That is, the air model of drawing 8 becomes the same as that of the model of drawing 7 , when the blow 
return from a cylinder has not arisen. 

In order to blow at the time of inlet-valve valve-opening initiation and for the local pressure buildup of the 
inlet port by the air current to also be reflected in an inlet- valve passage flow rate by accuracy by this, also 
in operational status with large bulb overlap, the restoration inspired air volume in a cylinder comes to be 
computed by accuracy. 
[0116] 

Drawing 9 is a flow chart which shows the outline of inlet- valve passage flow rate calculation actuation in 
which the inlet port model M25 of this operation gestalt was used. This actuation is performed as a routine 
repeated by ECU80 at intervals of [ comparatively short ] about 1ms. 

In actuation of drawing 9 , when it judges whether the blow return by the inlet port has arisen and blow 
return has arisen, the inlet port pressure which considered the local pressure buildup of the inlet port to an 
inlet port depended for blowing is computed, and an inlet-valve passage flow rate is computed using 26 
inlet-valve model formulas based on this inlet port pressure Pp. Moreover, when blow return has not arisen, 
an inlet-valve passage inhalation-of-air flow rate is computed, using several 25 as Pp=Pm and Tp=Tm. In 
this case, several 25 becomes the same formula as several 20. 
[0117] 

That is, in drawing 9 , the inlet- valve passage inspired air volume mc is computed using the inlet port 
pressure Pp computed last time at the time of actuation activation at step 901 . 

And at step 903, it judges whether the blow return by the inlet port from a current cylinder has arisen. 
Although this judgment may be judged based on whether the inlet port pressure Pp computed last time is 
lower than Pc, in having judged the existence of return [ blow ] only by the comparison of a pressure, the 
error under the effect of pressure pulsation of an inlet port may arise. For this reason, when a minute of the 
inlet-valve passage flow rate mc is integrated several times as a negative value and this integrated value 
turns into a negative value in forward and the inlet- valve passage flow rate in Pp<Pc (several 26), he is 
trying to judge the inlet-valve passage flow rate in Pp>=Pc (several 25) that the blow return by the inlet port 
from a cylinder has occurred with this operation gestalt. 
[0118] 

At step 903, when blow return has occurred (i.e., when the integrated value of the inlet-valve passage 
inspired air volume mc computed at step 901 is a negative value), based on the inlet port model formula 
progressed and mentioned later, the inlet port pressure Pp is computed to step 907. 
Moreover, when the integrated value of the inlet-valve passage flow rate mc is forward at step 903 (i.e., 
when blow return has not arisen), the pressure-of-induction-pipe force Pm computed as an inlet port 
pressure Pp in the above-mentioned inlet-pipe model formula (from several 1 6 to several 1 8) is set up, and 
this actuation is ended. 
[0119] 

Since an inlet-valve passage flow rate is computed using the pressure-of-induction-pipe force, using the inlet 
port pressure computed by this based on the inlet port model formula when blow return had arisen at the 
time of next actuation activation when blow return has not arisen, the exact inlet- valve passage flow rate in 
consideration of the effect of blow return is computed. 

Therefore, the amount Mc of the new mind with which it fills up in a cylinder comes to be computed by 
accuracy by integrating an inlet- valve passage flow rate over a 1 cycle intake stroke in consideration of 
positive/negative. 
[0120] 

Next, the inlet port model used for inlet port pressure Pp calculation of the above-mentioned step 907 is 
explained. 

The inlet port model which can be used with this operation gestalt is not necessarily restricted to one, and 
some its inlet port models are usable. Although this operation gestalt explains the example of an usable inlet 
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port model below, in actuation of drawing 9 , an inlet port pressure is computed based on either among the 

following inlet port models. 

[0121] 

(1) The 1st operation gestalt 

Drawing 10 is drawing showing the inlet port model of this operation gestalt typically. 
With this operation gestalt, when blow return has arisen, an inlet port is assumed to be the closed space 
which has the predetermined virtual volume, and it blows over this closed space from a cylinder, and asks 
for the pressure buildup at the time of presupposing that the gaseous whole quantity flows, and the pressure 
which added a part for this pressure buildup to the pressure-of-induction-pipe force is used as an inlet port 
pressure Pp. 

blowing — if a high-speed air current flows into a short time immediately after inlet- valve valve opening 
from a cylinder in an inlet port, in an inlet port, lifting of a pressure will arise locally with the flowing gas. It 
is thought that this local pressure buildup becomes so large that a flow rate is so large that the flowing gas is 
high-speed. 
[0122] 

So, with this operation gestalt, an inlet port assumes that it is the closed space (container) which has the 
volume of the imagination which becomes settled according to engine operational status, and assumes that 
the pressure buildup which flows into this closed space and which is blown and is produced according to an 
air current is equal to the pressure buildup of an inlet port. 

Now, per **** and unit time amount blows the volume as a closed space of an inlet port, and pressure 
variation dPp/dt is expressed with the following several 27 from an equation of state per flow rate, i.e., the 
unit time amount of the inlet port Pp according [ will blow if it blows and the inlet- valve passage flow rate 
at the time is set to mc and ] to a gaseous inflow. 
[0123] 

[Equation 27] 



dPp 
dt 



Vp= mc RT 



[0124] 

Here, R blows, the gas constant of gas and Tp are the gas temperature in an inlet port, and it is thought that 
they are equal to the gas temperature Tc in an abbreviation cylinder. 

Therefore, inlet port pressure-buildup width-of-face deltaPp depended for blowing blows, becomes equal to 
the integral value of several 27 after initiation, and when change of Tc is disregarded, it is expressed with 
the following several 28. 
[0125] 

[Equation 28] 



RTc r 
APp = mc-dt 

Vp J 



[0126] 

Here, it asks by integrating the flow rate mc (blowing flow rate) computed by several 26, and the integral 

value of the inlet-valve passage flow rate mc is the inlet port pressure Pp, 

Pp=Pm+delta Pp 

It is computed by carrying out. 
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[0127] 

At the drawing 9 step 907, the actuation which computes the inlet port pressure Pp by the above, and the 
actuation ( drawing 9 , step 901) which computes the next inlet- valve passage flow rate mc using the value 
of computed Pp are changed serially, and the inlet- valve passage inspired air volume mc comes to be 
computed by ****** at accuracy. 

In addition, the virtual volume **** of an inlet port is beforehand computed based on an experimental 
result, for example, is stored in ROM82 of ECU80 as a numerical map using the inhalation-of-air valve 
timing VT of the engine engine speed NE and an engine 1, and reading appearance is carried out by ECU80 
based on a actual engine engine speed and valve timing for every actuation activation of drawing 9 . 
[0128] 

For example, it blows, so that an engine's rotational frequency becomes high, and a gaseous amount 
decreases. For this reason, the virtual volume **** is set as such a small value that a rotational frequency 
increases. Moreover, although it blows (so that bulb overlap is large) and a gaseous amount increases so that 
valve timing carries out a tooth lead angle, the pressure buildup of an inlet port blows and does not increase 
corresponding to a gaseous amount. 

For this reason, for this reason, the virtual volume **** is set as such a big value that the valve timing of an 

inlet valve carries out a tooth lead angle. 

[0129] 

(2) The 2nd operation gestalt 

The virtual volume **** of an inlet port was made into the constant value which becomes settled according 
to engine operational status (for example, an engine engine speed and valve timing) with the operation 
gestalt of the above 1 st. 

However, by actual operation, the inlet port is connected to the surge tank instead of a closed space, it did 
not flow into an inlet port, and does not blow and the whole quantity of a flow rate does not necessarily 
remain in an inlet port. Compressing the gas which exists in an inlet port or an inlet pipe, the gas which 
flowed serves as a pressure wave and advances the inside of an inlet port toward a surge tank. 
[0130] 

So, with this operation gestalt, the virtual volume of an inlet port assumes that it expands according to an 
internal pressure. Moreover, as shown in drawing 1 1 , in the case of expansion of the virtual volume of an 
inlet port, the resistance expressed with the spring S of spring constant k and the dash pot D of an 
attenuation coefficient c to a boundary with the existing gas in an inlet port arises, and it resists and blows 
over this resistance and it is assumed that a gas expands. 
[0131] 

Now, the law of conservation of energy of the gas which will flow into an inlet port if it blows, and a 
gaseous flow rate is set to mc and temperature is set to Tc which flows into an inlet port is expressed with 
the following several 29. 
[0132] 

[Equation 29] 



d 



^ 1 Pp-VpV — R-mc-Tc- dw 



dt^k-1 J k-1 dt 



[0133] 

Here, although W is work which flows into an inlet port and which is blown and a gas does, if expressed as 
work which pushes the movable common wall 1 1 0 of the imagination energized with Spring S and the dash 
pot D in this work as shown in drawing 1 1 R> 1, the following several 30 will be obtained. 
[0134] 

[Equation 30] 
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dw _ dx 

= F» — 

dt dt 



Pp*A-K 



Vp 



dt(Aj) 



dt 




[0135] 

Here, A is the port cross section (constant value). 

Several 3 1 will be obtained, if the several 30 right-hand side is developed and it places with k/A=K and 

c/A=C. 

[0136] 

[Equation 3 1 ] 



^=p P ^E-K.Vp.^E-cf d ^l 

dt v dt dt V dt ) 



[0137] 

Therefore, several 3 1 is substituted for several 29 and several 32 is obtained. 
[0138] 

[Equation 32] 



dPp k-R-Tc kPp dVp k-1 

• me — • — —+ 



dt Vp 



Vp dt Vp 



2\ 



dt { dt ) ) 



[0139] 

Moreover, the several 33 following relation is materialized between ****, Pp, and mc from an equation of 

state. 

[0140] 

[Equation 33] 



RTc r 
Vp= Imcdt 

Pp J 



[0141] 

Furthermore, if it is discretized and several 32 is expressed, the following several 34 will be obtained. 
[0142] 

[Equation 34] 
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[0143] 

The inlet port pressure Pp is computed by [ of several 34 ] calculating serially using the inlet-valve passage 
flow rate mc (blowing flow rate) computed by several 26, and several 33. The constants K and C in several 
34 are beforehand called for by experiment as the engine engine speed NE and a function with the inlet- 
valve valve timing VT, and are stored in ROM82 of ECU80 as a numerical map using NE and VT, and 
reading appearance is carried out from this map based on a actual engine engine speed and valve timing 
every several 34 operation. 

This is enabled to compute the inlet port pressure buildup depended for blowing over accuracy further, and 

the calculation precision of the restoration inspired air volume in a cylinder improves. 

[0144] 

(3) The 3rd operation gestalt 

Although the inlet port internal pressure Pp was computed with the 2nd above-mentioned operation gestalt 
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based on the law of conservation of energy of the work which the movable common wall 1 1 0 of imagination 
of drawing 1 1 does, and the gas which flows into an inlet port, in order to simplify count, the balance of the 
force which acts on the both sides of the movable common wall 1 1 0 of drawing 1 1 is considered with this 
operation gestalt. 



Since the pressure-of-induction-pipe force Pm and the force F by Spring S and the dash pot D are acting on 
the left-hand side (inlet-pipe side) of the movable wall 1 1 0 of drawing 1 1 and it is acting on the right-hand 
side (inlet port side) of the movable wall 1 1 0 in the inlet port pressure Pp, it will become deltaPpxA=F if it 
places with deltaPp=Pp-Pm. It is the sum total of the energization force according [ F ] to Spring S and a 
dash pot D according [ A ] to the area (cross section of an inlet port) of the movable wall 110 here, and is 
given by the following several 35 using spring constant k of Spring S, and the attenuation coefficient c of a 
dash pot D. 
[0146] 

[Equation 35] 



Using several 35 and above-mentioned deltaPpxA=F, if it places with K=k/A2 and C=c/A2, several 36 will 

be obtained. 

[0148] 

[Equation 36] 



= K-Vp+C-^ 

dt 



[0149] 

Moreover, **** and dVp/dt are expressed with a-37 number and several 38, respectively. 
[0150] 

[Equation 37] 



[0145] 




[0147] 




[0151] 

[Equation 38] 
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dVp RTc 

= *mc 

dt Pp 



[0152] 

deltaPp blows, and since it is the local pressure-buildup width of face of the inlet port by gaseous inflow, it 
serves as the inlet port pressure Pp. 

Therefore, first, deltaPp sets with Pp=Pm and ****=o as initial value at the time of inlet- valve valve 
opening, is calculated from several 38 using the inlet-valve passage flow rate mc (blowing flow rate) 
computed by several 25, substitutes this for several 36, calculates deltaPp, and calculates Pp from 
Pp=Pm+delta Pp further. And mc can be calculated by several 25 using the value of calculated Pp, several 
38 to several 37, ****, and dVp/dt can be computed using this mc and its integrated value, and it can 
compute by substituting for several 36 by serial count which repeats actuation of calculating deltaPp. 
[0153] 

(4) The 4th operation gestalt 

Next, the 4th operation gestalt of an inlet port model is explained. 

With this operation gestalt, it is assumed that the inlet port pressure Pp when blow return occurs becomes 
settled by the pressure loss at the time of duct passage of the imagination which blows and contains the inlet 
valve of an air current. 
[0154] 

Here, when the gas flow rate at the time of duct passage is set to u,.a pressure loss (Pc-Pp) serves as a 

function of u2, and is expressed with the following several 39. ^ 

[0155] 

[Equation 39] 



(Pc-Pp)=lpu 2 +f(u 2 ) 



[0156] 

Here, rho of the first term of the several 39 right-hand sides blows, is a gaseous consistency and becomes 

the cylinder internal pressure Pc and a function with temperature Tc. 

[0157] 

Moreover, the second term of the several 39 right-hand sides is expressed like several 40. 
[0158] 

[Equation 40] 



2 1 « 

f (u )=cf — -u 

d 



[0159] 

It is with a constant, and the die length of the passage of imagination and the diameter in which alpha 
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contains 1 in and d contains an inlet valve here, respectively, and becomes the function of the lift L of an 
inlet valve. About the value of alpha, 1, and d, it has asked by the experiment using a actual engine 
beforehand, and has stored in ROM82 of ECU80. 
[0160] 

Moreover, Pc and Tc use the value computed with the cylinder model ( drawing 7 , M40). Thereby, the inlet 
port pressure Pp is computed by the following several 41. In addition, inlet port pressure-buildup part 
deltaPp depended for blowing from several 39 can also be expressed as deltaPp=Pp-Pm=Pc-Pm-(rhou2/2)- 
alpha(l/d) u2. 
[0161] 

[Equation 41] 



Pp = Pc--p u 2 -f (u 2 ) 
2 



[0162] 

Namely, mc computed with this operation gestalt using the cylinder model ( drawing 7 , M40) and the inlet- 
valve model ( drawing 7 , M30), By repeating the actuation which computes rho and u using Pc and Tc, 
computes the inlet port pressure Pp using several 41, and computes mc, following Pc, and following Tc, 
local lifting of the inlet port pressure depended for blowing is reflected in the inlet- valve passage flow rate 
mc, and the calculation precision of the restoration inspired air volume in a cylinder improves. 
[0163] 

(5) The 5th operation gestalt 

Next, the 5th operation gestalt of an inlet port model is explained. 

Drawing 12 shows the inlet port model of this operation gestalt typically. As shown in drawing 12 , with this 
operation gestalt, it is assumed that it flowed into the inlet port, blow from an inlet valve, gaseous [ a part 
of] extracts, 120 is passed, and it flows into an inlet pipe. 
[0164] 

Now, a-42 number is obtained from a gaseous equation of state using the gas flow rate mp which passes the 
pressure-of-induction-pipe force Pm, the inlet port pressure Pp, the inlet port temperature Tp, the inlet port 
volume ****, the inlet-valve passage flow rate (blowing flow rate) mc, and drawing 120, and flows into an 
inlet pipe from an inlet port. 
[0165] 

[Equation 42] 



dPp RTp / v 

— — = — — (mc-mpj 
dt Vp 



[0166] 

moreover, the flow rate mp which passes drawing 120 — Pm, Pp, and Tp — and it blows, extracts using the 
gaseous ratio of specific beat kappa, and is computed by 43 formulas of a passage flow rate, and several 44. 
[0167] 

[Equation 43] 
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mp=Cp* Ap 



Pp 



^Pp> 



[0168] 

[Equation 44] 




rpmV^ ' 




I'm 



\ 



k+l 



IPPJ 



PmJ 2 V-i 
Pp = ^k + lj 



k N 



— <f— 1 



k-1 



[0169] 

Several 43 Cp is the flow coefficient of a converging section, Ap is the opening area of the virtual drawing 
1 20, and it is set here according to an engine's operational status (the engine engine speed NE and inlet- 
valve valve timing VT). 
[0170] 

Here, it blows, and if the gas temperature Tp rewrites several 42 to a discrete format, assuming that it is 

equal to Tc whenever [ abbreviation cylinder bashful temperature ], it will obtain several 45. 

[0171] 

[Equation 45] 



RTc 



Pp (t + At)= Pp (t) + (mc - mp)- At 



[0172] 

With this operation gestalt, it blows, the inlet- valve passage flow rate mc is computed from an inlet- valve 
model by making Pp=Pm and mp=0 into initial value using several 26 at the time of generating, and the inlet 
port pressure Pp is computed from several 45 using this mc. Furthermore, from next time, the actuation 
which computes several 43 to mp for mc from several 26 using this computed Pp, respectively, and 
computes the inlet port pressure Pp from several 45 using these is repeated. Thereby, local lifting of the 
pressure to an inlet port depended for blowing is reflected in the inlet port pressure Pp by accuracy, and the 
calculation precision of the restoration inspired air volume in a cylinder improves. 
[0173] 
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(6) The 6th operation gestalt 

In this operation gestalt, the inlet port pressure buildup after inlet- valve valve opening is approximated with 
a sine function. In actual operation, the inlet port pressure buildup after inlet-valve valve opening is not 
necessarily a sine wave, though natural. However, since it is next very short time amount at the time of inlet- 
valve valve-opening initiation that blow return becomes a problem, gross errors are not produced even if it 
approximates a pressure buildup with a sine function. 
[0174] 

This operation gestalt expresses pressure-buildup deltaPp=Pp-Pm after inlet- valve valve opening by several 

46 as a sine function of crank angle CA. 

[0175] 

[Equation 46] 



APp = 



asin(b-CA-CAio) CA-CAio ^ 7c 
0 CA-CAio> n 



[0176] 

Here, CAio is a crank angle to which an inlet valve starts valve opening. 
[0177] 

Since it blows and is dependent on a gaseous flow rate, both the amplitude a (maximum-pressure lifting 
width of face) and periods b (between pressure rise periods) of the sine function in several 46 become a 
function with the ratio (Pm/Pe) of the amount OL of bulb overlap, and an intake pressure Pm and an exhaust 
pressure Pe, and can be expressed like several 47 and several 48. 
[0178] 

[Equation 47] 




C ( (OL + C 




[0179] 

[Equation 48] 



b=f 2 ^,OL^ ^ C 3 (OL + C 4 ) 




[0180] 

Here, CI to C4 is a constant, and is beforehand defined by the experiment using a actual engine. 
[0181] 

Drawing 13 is drawing showing an example of inlet port pressure-buildup width-of-face deltaPp computed 
using several 46 to 48. In drawing 2 , curves VT40, VT20, and VT0 show the case where inlet- valve valve- 
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opening crank angles (several 46 CAio) are 40 degrees, 20 degrees, and 0 times in front of a top dead center 
(BTDC), respectively. 

In addition, several 46 the amplitude a and period b of a function may store OL (Pm/Pe) in ROM82 of 
ECU80 beforehand in the form of the numerical map made into an argument instead of expressing using a 
formula like several 47 and several 48. 
[0182] 

It is possible by expressing with this operation gestalt to several 46 form, assuming inlet port pressure- 
buildup width-of-face deltaPp after inlet-valve valve opening to be a sine function to become simply 
possible to blow over accuracy and to compute the inlet port pressure buildup at the time of generating, and 
to raise the calculation precision of the restoration inspired air volume in a cylinder. 
[0183] 

[Effect of the Invention] 

When computing the restoration inspired air volume in a cylinder using the inlet- valve passage inhalation- 
of-air flow rate which computed and computed the inlet-valve passage inhalation-of-air flow rate based on 
the inlet- valve model formula according to invention given in each claim, to accuracy, it can blow, a flow 
rate can be computed and the common effectiveness it is ineffective to the thing from a cylinder to an inlet 
port for which the calculation precision of the restoration inspired air volume in a cylinder is raised being 
possible does so by taking the pressure of an inlet port into consideration. 
[Brief Description of the Drawings] 

[Drawing 1] It is drawing showing the outline configuration of the operation gestalt at the time of applying 

this invention to a jump-spark- ignition type engine. 

[Drawing 2] It is drawing explaining a throttle- valve model. 

[Drawing 3] It is drawing explaining an inlet-pipe model. 

[Drawing 4] It is drawing explaining an inlet- valve model. 

[Drawing 5] It is drawing explaining a cylinder model. 

[Drawing 6] It is drawing explaining an exhaust valve model. 

[Drawing 7] It is the block diagram showing the whole air model in the conventional technique. 
[Drawing 8] It is the block diagram showing the whole air model to have an inlet port model. 
[Drawing 9] It is a flow chart explaining the inlet- valve passage flow rate calculation actuation in the air 
model of drawing 8 

[Translation done.] 
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1 .This document has been translated by computer. So the translation may not reflect the original precisely, 

2 **** s hows the word which can not be translated. 
3. In the drawings, any words are not translated. 



DRAWINGS 



[Drawing 1] 




[Drawing 2] 
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[Drawing 3] 



ml 
Ta 



Vm , Pm, Tm 



mc 
Tm 



R 



Pni<k>-Pm(k-1)-*- At* ^ (mt(k-l)Ta-ncA V E (k-)) Toi(k-O) 



Tn(k) 



Pra(k) 



[Drawing 4] 



http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_web_cgi_ejje 



6/8/2006 



JP,2004-278359,A [DRAWINGS] 



Page 3 of 7 



E 

Ql 




+ 

I 

i 



csj 



IS 
I* 



+ 



rr 



CL 



> 



< 



o 
11 

o 
e 



-- * 

y 



[Drawing 5] 
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[Drawing 6] 
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[Drawing 7] 




[Drawing 8] 
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[Drawing 10] 
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ftTei:*tttfii-*Sf$i;aft-t>-9-tSf^CE-feV-9-t**Ktt'bti.TUl3 ( II -f ) , * 

*if ftflfJtlKP e i**EjTe ki:*t£.t4f **ECU8 0O>f>?-7i-^8 

5 \ZXti Ltv>4, 

[ 0 0 4 5 ] 

^illltli, ECU8 0 ti H Kl 1 0Ov>; > *r z 0 i:5E«?ii4f^f Mc £ ffl v> T 

<?- <r> v 'j ^^tciJJtsiRfl-Hltif c £ , f c = K I xmc t L t IT ta -f z> o k I tifiira 

[ 0 0 4 6 ] 

S , ECUS Oliy'j >72 0i:^«$n4f«i < tt l*i :fc « ift * * ) Mc^^raio 

^t^^t4 0 £ ^ x Wfc L fc ^ T ^ ^ ;V ft # j£ £ ffl v> T || fcB "T £ » 
[ 0 0 4 7 ] 

*iiii.Tiis$jf)i:, *m\&m xtftfrm 2 o o i -3 i 6 3 5 ot i:tstLtx7 
^;aMc#ai^ss-ife^-f^fci6ic % s -r # & 2 0 0 1 - 3 1 6 3 5 o-*»cij»t*«rt£ 

[ 0 0 4 8 ] 

■T >W t? S L T , W m ft Z %L W) -f Z> ty. Hi <T> & Jj y iS ft , aiia^fttiitCJ: t)t»S|| 
K$t«l*I(-?E»£*L&i&fi;g: ( firj 1*3 ?E * i& ^ fi ) Mc SJtHltii^KLT^S, 
[ 0 0 4 9 ] 

( 1 ) n h ;b # x* ;W 

xn 7 h;uf€f-;HJ, « $K I: l±0 2 i:f t i t I: t 9 (*nyh*#BP») £ il i§ 
-f&^#:<Oj?fi*ti: LT4x.f, jx4„ 
[ 0 0 5 0 ] 
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* 7 )l> ff # T' ^ $ it £ o 
[ 0 0 5 1 ] 

[ft i ] 



[ 0 0 5 2 ] 
[ft 2 ] 



mt = Ct(9t) At(9t)P»JIi±rrjf -JE_V -(-?*_ -2_Y 

V2KRTm VU + lJ Urn K + lJ 



[ 0 0 5 3 ] 

£ £ T , Rlli^n? U#4 3ia« ( * « ) O JE 7J T s£ fill ( !ft *t * IB ) OjETjJ: 
•)Sv'%^(Oi1?J)K^S!;il (JHsfE) 73fo]<0^.n-7h^#iiii>S.^BftS : $r, ft 2 (4 * n 
h ^ # 4 3 J:a«OE*#Til[|ll©EJJ: »)fi<aa#4tTi-'4»#OXD7 h A- # il 
$ ^ 51 sit S £ > f ^?tLgLtv>4o 
[ 0 0 5 4 ] 

ft 1 , ft 2 l: £> v> T , Ct ( 0 t ) l± s!E * i& ft , At ( 0 t ) tt^n? l-*#|PIBS|-e 

* 0 , Ct ( 6 t ) , At ( 0 t ) li-fn?tl^0 7 >^#SSM(OHSfab5o Pa 
(ixn-y h*#±««*>BE:&-e**Jfc^Bl-C(i**CflE* N P m l± * n y h^#T»fl«E 

SR-CJiift*CtfE27fcSEl- 0 * *: , T a , Tmli, f it-ftixn 7 F^#±I 

■9 o 

[ 0 0 5 5 ] 

* * Jfc J£ Bt T 14 % iiiiC t ( 0 t ) S^^nyh^^RpBltAt ( 0 t ) 14 ^ *t> J* 
n ? h^#HKtf#i:ffa$it, &t'7'y7'Of TECU 8 OOROM8 2 K f& fft $ ;h, T 

ft 1 . S2 i:f f i ^ i:, $IHIioxn'; hi-ftf ;ni, *:^EP a S^K^f E 
TJPm, TC^fiST a i tlif rtiatTm, SV^n 7 f- ;P # H j£ 0 t £ A 7J t L T 
, X n y h^#HiRHiltm t £ ffi ts £ ~T & * ^ * n# 5$; -? $> z> „ 
[ 0 0 5 6 ] 

( 2 ) m. ^ * * >v 

X. O'iSST m * » fib -t h o 
[ 0 0 5 7 ] 

X n y h;p#4 3^e>©^:#3 2 C14t-y ? >^ 4 4£-&tri&^;iF4 1 »#f 5r Vm 



(10) 



JP 2004-278359 A 2004.10.7 



£ £ <6 , aXlRStailim t , I^SSliTahiJo 

[ 0 0 5 8 ] 

[ 0 0 5 9 ] 

fi£ o X , K S f£ iij >J» f> T 12 ft 3 rt s , i^^-Sff!l«'tTEa4i« < ft -T ft, 3M> ft 

& O 

[ 0 0 6 0 ] 

[ft 3 ] 



d (Pm / Tm ) / dt = (R / Vm) • (rot - mc) 



[ 0 0 6 1 ] 
[ft 4 ] 



dPm/dt = ic-(R/Vm) (mt -Ta-mc Tm) 



[ 0 0 6 2 ] 

ECU8 OtOfliiiiff^A, ft 3 , S4 * L Tlftft L it 

, T O ft 5 „ ft 6 £ ffl v> & „ 

[ 0 0 6 3 ] 

[ft 5 ] 



Pm P in 1^. 

— (k)= — (k- 1)+ At— (mt(k - l)-mcAVE(k - 1)) 
l m l m vm 



[ 0 0 6 4 ] 
[ft 6 J 



Pm (k)= Pm (k - l)+Ahc-5-(mt (k - l)Ta - mcAVE (k - l)Tm (k - 1)) 



[ 0 0 6 5 ] 

z z x , k a^mnn^ftftom* , k-i(±mriift#iifi ; p#?)<i£gL, a t i± 
ik it + & m <n & a n («xif8m S gto^ vva) jsltvSo t tz , ft 5 , ft6 

K tt & m c A V E (i : miEi-^ISiF'g^^^amc <D®.%fim. 1 * -f 
[ 0 0 6 6 ] 

^^iSJgTmli, S5 iD**tPm/Tmt, ft6^<b**^Pmi:*fflv»TTm = 
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P m / (Pm/Tm) tLTStBJiX^o 
( 3 ) R£ ^ # f ^ 

[ 0 0 6 7 ] 

&ot, m.M.ft* r n> amm<?) * v y \> *, # ^ ^ * (ai, » 2 ) ii«t:TE«t7, 

[ 0 0 6 8 ] 
It7] 



mc = Civ(L)Aiv(L)PmJ-^- 1 -JS-Y L_Y 

V2KRTm VU + lJ lPm K + lJ 



[ 0 0 6 9 ] 

ims] 



[ 0 0 7 0 ] 

$1 * , Pm, Tmt±ia^l : jE^S.O r iaK, Pc, T c 14 -> V > 7 ft £E 2j £ tT it £ h Z> „ 

* , Civ (L) l±»iES#S:T*l9, IStO'J 7 Aiv (L) 14 © £C 

#3 2?>HPiB!*-e*»K HUC^HLOBftUi. 
[ 0 0 7 1 ] 

* * ifc ff£ B T 14 , »f C i v (L) • A i v (L) Wttli^L WlC^LTtf»?it> Sfe 
if-y*OitECU8 OOROM8 2 Kfttt S JlT v» 4 „ 

x n y h * # * t* )u (ifei, sfc 2 ) t h tttc & 7 amM'g&.Jj p m a* -> »; > ^ 1*1 je # p c 

[ 0 0 7 2 ] 

Sfc7. S8i:ftiii:, ^IffiiSoS^^tf ^(i, 5i ^ EE 2j P m , v > ^ |*| EE 

Kitiat^itm c £ m £ + & o 
7>^^i8o° izttm-tzftmr (is oca) xm l tzmx $> z> 0 

[ 0 0 7 3 ] 

1~ & fc> *> , mcAVE=Mc/T (1 8 OCA) 
( 4 ) v > ^ 7* H> 

* £ ifc m <r> -> > y =e -r ;H4 , 1 5 c^ti 9 liftf v c ( *r ^ ) 
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[ 0 0 7 4 ] 

1- ft *> *> , *IiiS«y 'J >f*r*^fl±, m c , me, T m , T e ( R Xf * y > ? ft 
) £A7JtL, v 'J > fr ft S. 7J P c , ffi. $t T c #. If* ftfi 1 « + U 1*1 K ft * $ ft q£ ft o 
ij£. fi M c tffiAt So 
[ 0 0 7 5 ] 

v > ^ t-* ;n± , * * ft ic (i i * )u ^ & # m iz m o' < t IE ft 9 t m S ft * o 

[ 0 0 7 6 ] 
[ft 9 ] 



dPc tcRTm KRTe kPc dVc Ow 

— — = — — — "mc+ -me • + (k - 1)— — 

dt Vc Vc Vc dt V ; Vc 



[ 0 0 7 7 ] 

[ 0 0 7 8 ] 

* jfc , y'j^fliaSTcli, % # <D # 81 ^ S 5$; ^ <b T IB ft 1 0 T £ $ ft 4 0 

[ 0 0 7 9 ] 

[ft 1 0 ] 



Tc= (PcVc)/(MclR) 



[ 0 0 8 0 ] 

£ £ -C , Mc <M#fM#B#^<b3i:£ (T eft 

»B#) * -TOift*#fc oaEAfcflsa* (mc+me) 0»&ItLT4UJi 

^) o 

[ 0 0 8 1 ] 

± te ft 9 (i , :©ii©«?ttECU8 0?»*iffi:t#t!lftv^i6, JK ft -fb U 

[ 0 0 8 2 ] 
[ft 1 1 ] 
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[ 0 0 8 3 ] 

£ tes ft 1 l z & t i*$l9 *mtkit-$- z>l£<?>&frmwi ( v >; > y * ?> )i, -c i± ® % <g 

^ ^ tf) J# U t£ T A t t± ffi < |2 5g £ ft , Mx.lflmsISO^ 

T * & o 

I 0 0 8 4 J 

i'- # ® ii eft fi m e (i , *i:R« t 4 f * C J: ») ft 4 „ 

(5) »f SC # * 7* ^ (07, M5 0) 

jv {±WL%ft* t )v tmm^, » at # - v <ds.1j tsjs, p e , t e t>* -> >; > 

[ 0 0 8 5 ] 

II 6 (± , m% # * r-' )\> £ Tf; f » SI 6 tc ^ ■*- <fc ^ tc ^ ^ # * 7" ;Hi g * ift tc (± © ^ # r' )V 
tm-X-h*) , #|^#«i§5iL*l±ftl 2 $ ft <± ft 1 3T-4x.fjft4o 
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[ 0 0 8 6 ] 

i ft i 2 ] 



me = Cev (L) • Aev (L) • Pe. 



K + l 

2icRTe 




^K + l 




[ 0 0 8 7 ] 
[ft 1 3 ] 



me = Cev (L) • Aev (L) • Pc. 



K + l 

2kRTc 



/ Pe 



K + l 



1 

P c 
ft Z> 
e v 
L <D 
[ 0 

ft J3 

-t y 

tz tb 

x m 

i tz 

i> m 

I o 



2 (i „ Jft $i - h tf> <b v 

h EE7J P e «t 
# ( P c tf* P e £ *) 



IB] 



( L) l±*ft-F*L#fjJC# 

n ft t ft * o 

0 8 8 ] 



v > y ftiz m # & -t & 

•) <£v>»£-) , ft 1 3 14 v 
v»»^-) * -f- ti -T ft. 3? L T 

3 5 <d »iifi#ft t fflnmm 



& ( -f ft t> *> , ->'J >yfiEJ 
V > r <b #E 5t # - h K #f St tS j % 
v» * o , Cev (L) , A 



•f* ^ 

SC it 
[ 0 



ft inn. Jj£ 

0 8 9 
(4, _h 

M c <r> 

0 9 0 



ife B M x 

9 i <h *> 

p » a 

■fc > -9- £ 

] 

7* A' H * 

] 



14, JftjR; 
> "t ( HI 

xmfts. 

°T 1ft T* h 
ST e li 



ax 



L tz 7. D 

a* <b «* 

O T* ^ 

f C o v> 



# - bff^p e t m ft m. m 

i tit Hi -1- & ft t> *) \z P e * 

•y h A/ # -=E -f ;p „ CR $t @ 
tf**, A# k its 2k 



T e t4 m ft # - b tc K It tz ffl ft EE 
L^it*, t^lf T e i:ov>t 

t-* . is §t ^ ^ t* a- , -> >; > ^ 

■f^uyf&mx&Zo i 6 Irffl 



( A ) X n y V )V ft J* )V (07, Ml 0 ) 

( i ) nn 

[ 0 0 9 1 ] 
[ft 1 4 ] 



mt = Ct(9t)At(et)Pa. 



K + l 




Y 



2icRTa VVK + lj 




[ 0 0 9 2 ] 
[ft 1 5 ] 
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mt = Ct(0t)At(et)Pm ' K + 1 



2kRThi 




I 0 0 9 3 J 

( 2 ) A7J 

■ ^n? h*#MS j t 

• *SEP a (*SlE-tV*6 3 fflj) 
•±SCi&j£Ta (BR^fflK -t > -9" 6 2 tiS 7J ) 

• t^f P m 

•t«ffitTm (S5cl 5 ( i& »fE B# ) <7)®&) 

( 3 ) 

• * n y b^#aaiR$taE*m t 

(B) l& % m * 7* ^ (17, M2 0) 
( 1 ) 1 1 * 5t 
1 0 0 9 4 ] 
[ft 1 6 ] 



10 



20 



—00= — (k- 1) + At — (mt (k - l)-mcAVE(k - 1)) 
Tm Tm Vm 



[ 0 0 9 5 ] 
[St 1 7 ] 



30 



Pm (k)= Pm (k - 1)+ AtK — (mt (k - l)Ta - mcAVE(k - l)Tm (k - 1)) 



[ 0 0 9 6 ] 
[ft 1 8 ] 



Tm(k)= Pm ^ pin 

fm" (k) 



40 



[ 0 0 9 7 ] 

( 2 ) XtJ 

• xu -y \ )V i?mi&W.%.ffi.mm t (k - 1) 

• (¥-*HVL) mcAVE (k-1) 

•Pm(k-l),Tm(k-l), ( P m / T m ) (k-1) 50 
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( 3 ) &D 

■ W.^.WI±t) P m ( k ) 

• w.m.WM.m t m ( 

( C ) ««#tf^ 
( 1 ) ft" » 5* 
[ 0 0 9 8 ] 

im i 9 ] 



k ) 

(1 



M 3 0 ) 



mc = Civ (L) Aiv (L) Pm 



K + l 



2KRTm 




K + l 



-{—-—) 

{Pm K + l) 



10 



[ 0 0 9 9 ] 

[Ifc 2 0 ] 



mc = Civ(L)Aiv(L)Pc 



K + l 

2kRTc 




20 



[0100] 

( 2 ) XJl 

t=JE tl P m 
®.%Wi&m T m 
-> > V J± tS P c 

>^rtSST c (12 0 (iS»£B#) <0 
^ * 7* 'J 7 ML ( ? 9 > ^ £ h M m ) 

(3) aiTj 

W. % ft m & W. Hi. %L M (BftRftt) mc 
• KI#IiifMltmc AVE (¥*^ffi) 



30 



( d ) ~s >) > v ^ -r ?v 
( i ) ftl?^ 

[0101] 
2 1 ] 



(07, M 4 0 ) 
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[ 0 1 0 2] 
[&22] 



40 



Tc=(PcVc)/(MclR) 



I 0 1 0 3] 

( 2 ) XtS 

• ~> *) > y&mv c it ? > ? n& hnm) 

• i&^ifrit SgK^S (Rtttt) mc 
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( 3 ) mti 

■ > 'J > ^ EE 7J P c 

■ '/T/^rtitTc 
• fgi 1*1 5E * S M c 

(E) $ 41 # * f * (07, M5 0) 
( 1 ) It » A 
[ 0 1 0 4] 
[t23] 



me = Cev (L)- Aev (L) • PeJ K + 1 . If-^-l - f — — 1 * 

V2kRTc VU + lJ lPe K + lJ 



[0105] 
1824] 



me = Cev(L)-Aev(L)-Pc 1 /-^- Jf-^-T -f — — V 

V2icRTcVU + lJ Uc K + lJ 



[0106] 

( 2 ) A7J 

• Sf^fEAP e 

• HIKflKT e 

• -> V ^ EE tS P c 

• -> 'J >^rtiSST c (ft 2 4 «> 4ft g- ) 

• /< ;l/ r 'J 7 ML (?7>?£§^t>#ffl) 
( 3 ) ffl 7J 

J#»#fl«»&SlE* (Rfiftt) me 

± & O X. -9 , * n y Y )U # r-* ;l/ t p # •=£ 7* ;w JSL O ^ , 1" & *> *> iR & § ^ -f )V 
[0107] 

K {C , ^nyh^#*^^OA*4, XV?h)Vftmm.8 (t) , ^tffiP a, * :SC iS jg 
T a £ E£ § , iotf ^0(BJtt#ttffl"P*4, fit, *IffiP a, ^^iStTalil 

«£oT, 17 i:f LfcxTtfii-f m$fi(ii:xD h frfifflm e (t) £ IS* Kl I! 4E ft N 

o 

[0 1 0 8] 

*#»Ctt, «Atf«Hft»*i:iia)miiLTPn, Pct^IEPaC, T m . Tc 
S'^KiSST a I:, £ x cSfimt, mc (±0 K-f-Jv-fftEifeL.. XD 7 ( 

t) Rv ! ?7>?nzm^x±&(7>mi4frb2 4<7)nm:fti (p e> t e « n a a * 

+ 4 ) o 
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^ LT, 3 <D $g 3M# h tL fz P m , PcTm, T c , m t x mc, mef OffffiUTil 

4 *> h $k 2 4<7>*EoU-2r&fMo *LT, z 0) ft ft i & M ft & ft & e> ft & & ft ft ( #fl ^ 

lf8mslS) -C -5" <T) B# jft U j3 It -S 0 (t) , NEc0ttt|trIlC0ft-»T*5fci6 7t:Pm, t m 
, mt, m c , m e f <?)f ^ffl^tl *) o Z ix IZ X *) , p m> tm, mt, mc^li 

ft»ga^f*SB#^r-^|^offl(;jRmt, UROflffrtftiJliRStic: - «t K & £ (i -r 

■C * £ o 

[0109] 

±EWg*<7)x7tf JVTM^ f^#^f^ ( EI 4 ) i$ V» T , v 'J > r t> IJft ft sK - h 
^<r>&.ffi.-t> t $L£Z,W,i? ( » 2 0 ) £ & # lit L £ ® ft 4 ft 5e ^ ft T v> & , 

*7t-/<5 ? 7 * * § < T^*w?Jc§jiL3j«^c* i o fti*i§i-lc&fi5i*!3fciil'Rft»ti 
jE5SUJt5fe$tL4 lit TJbi, 

ft * O * 5c ffl St ?> f T f & M a tJ* £. U £ o 
[01 10] 

* 36 W V 14 , lKI*Oi7tf JH:itLt, ®ftt=fcq&fttf-bi£-55-ftT*7'>MI:L 
, i& ft *° - F^f^i-fflv^ti^liLlltlilwt^-f- ^7j£i&ftt : JE7Ji:liSiJlI<7>« 

0 **»§9,<0x7*7*O^ftfiSS:|ill8t4H7 tHI47D ^lf*i 0 

[01 1 1 ] 

@ 8 c^t <t ? ^ m m m m t- f± ft « * r >i> m 2 0 t m. ft # * r * m 3 0 t <n n t- , 

ift ft # - Kf^M2 5 IS It <b ft T v» A £ ffi jg 4 „ 

R»# - h€7*M2 SOffiCovTttftatS*', # fl flj S "C t± .1 O IR ft * - b ■=£ 
r )VM 2 5 & ffl v» T m. ft # - hJETJP p £ g tH L , ®^#iii§gftftmc £ tfc 2 0 <7> ft t> *) 
C #C K 5? -f & 2 5 , »2 6 *fflv>Tl(iit4 J: a i:LTv>4, I ; f , Sfc 2 5 li Si ft # - 
b-^<7)PK§jgL75 s £i;Tv>&v>Ji§# ( -f *> *> , Pp^Pcco^-g-) , gt 2 6 f± 95. ft - 

[0112] 
[Sfe 2 5 ] 



mc = Civ(L)Aiv(L)Pp,pi±i- Jf _f * _ 

V 2KRTm VU + U lPp K + lJ 



[01 13] 
[St 2 6 ] 



mc = Civ(L)Aiv(L)Pc,/^±r f _JL_Y_f fE__J_Y 

V2kRTcVU + iJ lPc K + lJ 



[01 14] 

£ £ -C , f£ & & i 9 I- * ^ O q& ft - h ^ 7* ^ M 2 5 f IJ v 'J > / i> t> O P?: i i 

L#4CTV>4i^t:ii ( ft W JE 7J P m t (4 SUH K % P m «t A v» JE7J t L T W. ft # - 
h E 7J P p £ # ii} 1" & , # jf L C T v> & v» 1§ £• (4 , ft * - f EE # P p t LX 
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■f Z> Z t li ft t> & V> o 

[ 0 1 1 5 J 

7* )V t m — t( Z, o 

^#iijSsfEfiUiE5S(-S?!li$tt* J: -9 iz % & fz #> ^ /<? ;i> 7* * - 7 y -f -fr § v> a j£ Is 
»ciJv»-c<>jEauttrt^l<l%5C*35*J([ffi§ii* i ■? tc & & 0 
[0116] . 

it70-ft-l-TU t * 8 ftHi , ECU8 Ot:J:0«Alflm8fitOJtttftfflt> 
HI*-C«*«3ig3ft*A'-*-> t L X fi t> ix Z> 0 

m 9 O JA ft T- li , ?R ^, - b ^ O vfr § ig L Sl C T v» 4 5 jj» £ flj jg L > ift § il L rt* £ C 
T V> & % K li „ tyL%#- V <F> Vk & L \Z £ Z, - Y <V MIR W) % St) ±.g- Z iNVfc L 

m. * # - h J± 7J £ # ft L , £ tf> R *t # - h EE 77 P p K * o* v> T IR # ^ 7* ^ ft # 5*i & 2 
6*ffiv*Tift2C#a«ij|E*«:«aii-&o tfc, RSilLi»ti:Tv>4K*^i:ii, Pp = 
Pm. Tp=Tm<!:LT, »2 5*Jflv»-Cift*C#aiaiR*C8t**JHiii-,& 0 Z O ±§ IC li 
S5: 2 5 li ifc 2 0 tRtSKSii. 
[0117] 

1" & fc> *> , 0 9 i: i5 t ^ f 7 7' 9 0 1 TMi iitr ® ft ft It fr 1$ I- g as L £ <S ft - h EE 7J P 

■€■ L T , X-r-^ 7*9 0 3 tliSSy'i > ^ ^ <b % £t # - h I 1 L #4 C t v»i 

^JflttJ, i (O ® li , !W0*ffiLfcR*;K-l»EE2jPp;&tpcJ:»)te<fcoTv»* 
^ 5 tJMC » o* v> T f!l Bf LT & £V'i { , #ICj±7JOJtKO«^.-CP^§iI L O * * £ fij »r L £: O 
T li IR & # - hOE7J)!Ri*(0f$WlCj;^^Md ? feU^^i^-* f *^ o iofcft, $gft^§ 
TMiPp^PcO«r^Oia^:#ai§5fEg (12 5) 1:1, Pp<PcO^£-<OlRft#a:i§ 

atfi (t2 6) *n<7>m t lt, 85c0^o^*:#aig»E*mc hil, z nmm&tfn 

[01 18] 

^T?y9 0 3 t, l^llL^f 4Ltv^S^, t^^*>Xf 7-/9 0 Iff IliLtK 
^#ajg©^;4mcC7)aff®^ft<7)ffiT-$)c.^*§^|Cli, X f 7 7* 9 0 7 C I A , f£ i£ i" 

* « *C # - h ^ T"* ^ ft 3t 5Uc » o* v> T IR $t - 1* EE 7J P p illiitS. 

* tz , 7 7*9 0 3 tf ^#filSfm c <?)f 1" ig 

L#4tTv>4v>»^i;ii, "R ft ** - h EE 7J P p b L X m & <F> * 7 )i< ft nit. ( Sk 1 
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